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Abstract: In the context of global warming, drought events occur frequently. In order to better 
understanding the process and mechanism of drought occurrence and evolution, scholars have dedicated 
much attention on drought propagation, mainly focusing on drought propagation time and propagation 
probability. However, there are relatively few studies on the sensitivities of drought propagation to seasons 
and drought levels. Therefore, we took the Heihe River Basin (HRB) of Northwest China as the case 
study area to quantify the propagation time and propagation probability from meteorological drought to 
agricultural drought during the period of 1981—2020, and subsequently explore their sensitivities to 
seasons (irrigation and non-irtigation seasons) and drought levels. The correlation coefficient method and 
Copula-based interval conditional probability model were employed to determine the drought propagation 
time and propagation probability. The results determined the average drought propagation time as 8 
months in the whole basin, which was reduced by 2 months (i.e, 6 months) on average during the 
irrigation season and prolonged by 2 months (i.e, 10 months) during the non-irrigation season. 
Propagation probability was sensitive to both seasons and drought levels, and the sensitivities had 
noticeable spatial differences in the whole basin. The propagation probability of agricultural drought at 
different levels generally increased with the meteorological drought levels for the upstream, midstream, 
and southern downstream regions of the HRB. Lesser agricultural droughts were more likely to be 
triggered during the irrigation season, while severer agricultural droughts were occurred mostly during the 
non-irrigation season. The research results are helpful to understand the characteristics of drought 
propagation and provide a scientific basis for the prevention and control of droughts. This study is of 
great significance for the rational planning of local water resources and maintaining good ecological 
environment in the HRB. 
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1 Introduction 


Drought has become one of the most common natural disasters worldwide (Mishra and Singh, 
2012; Liu et al., 2018), typically inducing negative effects on the economy, agriculture, society, 
and environment (Cheng et al., 2013; Gu et al., 2020; Dahal et al., 2021; Zhang et al., 2022a). The 
economic losses caused by droughts during 1987—2006 in China accounted for more than 50% of 
meteorological disasters (Wang et al., 2012). It was estimated that under extreme drought 
conditions, the probability of corn yield reduction in Shandong Province of North China was 
greater than 80% and that in Jiangsu Province of East China was even greater than 86% (Zhang et 
al., 2022a). With global warming, drought assessment has become increasingly important for the 
establishment of drought prevention and early warning systems. 

The American Meteorological Society classifies drought into four types: meteorological 
drought, agricultural drought, hydrological drought, and socioeconomic drought (Heim, 2002). 
Meteorological drought generally occurs earlier than hydrological drought and agricultural 
drought. Once meteorological drought happens, the other two drought types (hydrological drought 
and agricultural drought) can occur easily (Huang et al., 2015; Guo et al., 2020). Although there 
are differences in the occurrence time of droughts, the direct impact and disasters caused by 
droughts are often manifested through agricultural drought and hydrological drought (Guo et al., 
2020; Zhou et al., 2021b). The propagation of water deficit signals between different types of 
droughts is defined as drought propagation (Apurv et al., 2017). Currently, research on drought 
propagation mainly focuses on propagation time, propagation probability, etc. 

Drought propagation time is usually derived from the correlation analysis between two types of 
droughts (Abbas et al., 2021; Um et al., 2022). It is determined by the maximum correlation 
coefficient between a drought type on an n-month time scale and another drought type on a 
1-month time scale; generally, the accumulation period n is defined as the drought propagation 
time (Li et al., 2016; Bayer Altin and Altin, 2021; Liu et al., 2021b; Xu et al., 2021). Based on the 
correlation method, some scholars found that propagation time from meteorological drought to 
agricultural drought in the Yellow River Basin was 0-3 months (Wang, 2020), and propagation 
time from meteorological drought to hydrological drought in the Pearl River Basin was 2-6 
months (Zhou et al., 2021b), both in China. Significant spatial and temporal differences have also 
been determined in drought propagation time (Huang et al., 2017; Li et al., 2020; Xu et al., 2021; 
Zhou et al., 2021b). For example, Xu et al. (2021) evaluated propagation time from 
meteorological drought to agricultural drought and found that it gradually increased from 2 to 7 
months from South to North China. Huang et al. (2017) concluded that propagation time from 
meteorological drought to hydrological drought had noticeably seasonal characteristics in the 
Weihe River Basin, China, with a shorter time in spring and summer, and a longer time in autumn 
and winter; while Li et al. (2020) determined a longer duration in spring and winter, and a shorter 
duration in summer and autumn in the upper Shaying River Basin, China. Ding et al. (2021) 
revealed that in North China, there was a stronger propagation relationship from agricultural 
drought to hydrological drought in summer and autumn compared to spring. 

Drought propagation probability is generally determined by either conditional probability 
together with copula functions or by Bayesian network model (Shin et al., 2018; Sattar et al., 
2019; Xu et al., 2019; Sun, 2021). Shin et al. (2018) employed the former method and found that 
propagation probability from meteorological drought to hydrological drought in Korea was 33% 
and 48%, based on standardized precipitation evapotranspiration index (SPEI) at the 3- and 
6-month time scales, respectively. Sattar et al. (2019) used Bayesian network model to predict the 
propagation probability of hydrological drought triggered by meteorological drought in Korea and 
determined propagation probability to range from 27% to 60%. The propagation probability of 
different drought levels has attracted attention over the past few years. For example, Dehghani et 
al. (2019) found that extreme and severe meteorological drought led to hydrological drought with 
a 70% propagation probability in the Karun Basin of Iran, while other levels of meteorological 
drought resulted in a normal hydrological drought with less than a 50% propagation probability. 
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Jehanzaib and Kim (2020) employed Bayesian network model and Markov random fields to 
examine the propagation probability of meteorological drought and agricultural drought to 
hydrological drought in the Nakdong River Basin of Korea, and revealed that the joint 
propagation probability from severe and extreme meteorological drought and agricultural drought 
to severe and extreme hydrological drought ranged from 49% to 83% on average. Xu et al. (2021) 
concluded that the propagation probability of agricultural drought in China increased in parallel 
with the severity of meteorological drought. Zhu et al. (2021) detected the propagation probability 
of soil moisture drought to be on the rise in most regions across the globe with increasing 
meteorological drought levels. 

The aforementioned studies have made progress in detecting drought propagation time and 
propagation probability across different regions; meanwhile, several scholars have also focused 
on the seasonal sensitivity of drought propagation time and the drought level sensitivity of 
propagation probability (e.g., Huang et al., 2017; Dehghani et al., 2019; Jehanzaib and Kim, 
2020; Li et al., 2020; Ding et al., 2021; Xu et al., 2021; Zhu et al., 2021). However, further 
research is still required to enrich the study of drought propagation law, especially for drought 
propagation probability and its sensitivity to seasons and drought levels. 

The Heihe River Basin (HRB) in Northwest China is characterized by sparse and concentrated 
precipitation and intense solar radiation, making the local drought problem prominent. Both 
drought and its propagation across the basin not only aggravate the local fragile ecological 
environment but also affect the water security, food security, and socio-economic development. 
Taking the HRB as the study area, this research aims to explore the propagation characteristics of 
drought propagation from meteorological drought to agricultural drought from two perspectives, 
namely, drought propagation time and propagation probability, and analyzes their sensitivities to 
seasons and drought levels. This work acts as a scientific basis for the management and planning 
of water resources in the HRB, and provides decision-making support in the early warning and 
prediction of droughts in this area. 


2 Study area and data collection 


2.1 Study area 


The HRB, the second-largest inland river basin in Northwest China, is located in the arid zone of 
Northwest China (37°44'-42°40'N, 97°37'E-102°06’E). The basin has an area of 142,900 km?, 
and can be divided into the upstream, midstream, and downstream regions from the south to north 
according to the locations of Yingluoxia and Zhengyixia hydrological stations (the hydrological 
stations at the outlets of the upstream and midstream of the HRB, respectively). The upstream 
region of the basin is the runoff producing area, and 90% of the water in the basin comes from the 
upstream (Wang, 2019). The midstream region is a runoff utilization area and is one of the top 10 
commodity grain bases in China. The cultivated land area accounts for nearly 20% of the total 
area of the midstream region. It is also the essential ecological safety barrier and transportation 
hub along the Silk Road Economic Belt. The downstream region is a runoff consumption area 
with a vulnerable ecological environment. 

The multi-year average annual precipitation and air temperature over the basin is 122.0 mm and 
5.1°C, respectively, which were estimated from the observed data at gauge stations of the basin 
from 1981 to 2015 via the Thiessen polygon method. Precipitation gradually decreases from the 
upstream region to the downstream region, while air temperature increases from the upstream 
region to the downstream region. Due to climatic characteristics and water intake for irrigation in 
the midstream region, the downstream region usually suffers from droughts and water shortages 
(Feng and Su, 2020), which makes the ecological environment of Ejina Oasis, a national ecological 
reserve in China, more fragile. Topographically, the upstream region has high elevations, with 
peaks above 4000 m. The midstream region is relatively flat with elevations from 1200 to 2000 m, 
and the downstream region is open and flat with elevations below 1200 m (Fig. 1). 
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Fig. 1 Overview of the Heihe River Basin (HRB) and locations of hydrological stations. Black solid line is the 
basin dividing line of the upstream, midstream, and downstream. DEM, digital elevation model. 


2.2 Global land data assimilation system (GLDAS) 


GLDAS is a joint development of NASA's Gotha Space Flight Center and NOAA's National 
Center for Environmental Prediction, aiming to simulate land surface fluxes and storage of water 
and energy (Rodell et al., 2004). The GLDAS has been updated to version 2 (GLDAS-2), which 
includes 3-h, daily, and monthly products with spatial scales of 0.25° and 1.0° (https://earthdata. 
nasa.gov/). GLDAS-2 also comprises the versions of GLDAS-2.0 and GLDAS-?2.1, in which the 
former covers the period of 1948—2014 and the latter covers the period of 2000-2020. Since both 
GLDAS-2.0 and GLDAS-2.1 involve the dataset from 2000 to 2014, we took the variable of 
precipitation as an example and compared the overlapping portion of the data in the HRB (Fig. 2). 
It is clear that the two curves representing each version are essentially consistent, with only slight 
differences in the extremes of individual years. Thus, the two datasets are largely comparable and 
consistent. Therefore, the data from GLDAS-2.0 during the period of 1981—1999 and the data 
from GLDAS-2.1 during the period of 2000-2020 were used in our study. 
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Fig. 2 Comparison of precipitation data derived from GLDAS-2.0 and GLDAS-2.1 during the period of 
2000-2014. GLDAS, Global land data assimilation system. 
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There are inevitably some uncertainties involved in GLDAS datasets (e.g., Liu et al., 2021a; 
Wt et al., 2021; Zhang et al., 2021a; Zhou et al., 2021a). For example, Wu et al. (2021) found that 
GLDAS-2.1 slightly underestimated the climatology values in northern China, and slightly 
overestimated the values during summer and autumn in southern China. However, GLDAS-2 
provides us with a substantial amount of available research data, particularly in areas where 
measured data are scarce. Several studies verified the reliability of GLDAS-2 dataset in specific 
regions. For example, Liu et al. (2019) used the Noah products from 1982 to 2015 in GLDAS-?.1 
to investigate the spatial-temporal characteristics of the dry-wet regime and vegetation dynamic 
responses in the Yarlung Zangbo River Basin in China, and showed that the spatial-temporal 
patterns of precipitation and air temperature simulated by GLDAS-Noah fitted well with those of 
the in-situ data. Li et al. (2022) used the Noah products in GLDAS-2.1 to investigate the 
vegetation dynamic variations and their associated responses to environmental changes in the 
Yarlung Zangbo River Basin. Zhu et al. (2021) compared the GLDAS-2.0 data with observed data 
from global sites, and found that the data from GLDAS-2.0 were objective and reasonable for 
global drought propagation analysis. In order to assess the applicability of GLDAS-2 data in our 
study, we compared the precipitation and air temperature data derived from both GLDAS-2 and 
gauge stations of the HRB over the same period (1981-2015). Results show that GLDAS-2 data 
was essentially consistent with the observed data based on spatial trends; the performance of 
GLDAS-2 data was better on the regional average scale, although precipitation and air 
temperature in the upstream region were underestimated. Therefore, GLDAS-2.0 and GLDAS-2.1 
data were selected as the basic data for our analysis. 

More specifically, data of precipitation, air temperature, and soil moisture from 1981 to 2020 
provided by Noah land surface model (Noah-LSM) in GLDAS-2.0 and GLDAS-?2.1 
(https://earthdata.nasa.gov/) were selected as our data sources. Table 1 reports the key information 
of these data. Soil moisture data include values at the depths of 0—10, 10—40, and 40-100 cm. The 
data were pre-processed prior to the drought analysis, including re-projection, conversion of units, 
etc. Figure 3a—c presents the spatial distribution of precipitation, air temperature, and soil 
moisture on the multi-year average scale. We then calculated the drought index based on the 
processed data (Vicente-Serrano et al., 2010; Martinez-Fernandez et al., 2015). In addition, land 
use and land cover change data in 2015 were provided by Landsat 8 (https://landsat.gsfc.nasa.gov/ 
satellites/landsat-8/). According to the current land use classification (Ministry of Agriculture of 
the People's Republic of China, 2017), there are 23 types of land use and land cover in the HRB, 
including forest, low-coverage grassland, cultivated land, desert, Gobi, etc. The distribution of 
each land use and land cover type is shown in Figure 3d. 


Table 1 Key information of the data used in this study 


Data Unit Spatial and temporal scales Data source 
Precipitation kg/(m°-s°) 
7 0.25°x0.25° . 
Air temperature K Monthly https://earthdata.nasa.gov/ 
Soil moisture kg/m? 
Land use and land cover change / 30m https://landsat.gsfc.nasa.gov/satellites/landsat-8/ 


Note: "/" denotes dimensionless. 


3 Methods 


3.1 Framework 


This research aims to explore the propagation time and propagation probability of meteorological 
drought to agricultural drought over the HRB. Regarding drought propagation time, besides its 
seasonal sensitivity, the dynamic characteristics were also considered since propagation time may 
not remain stable owing to the effects of climate change and human activities (Liu et al., 2021b). 
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Fig. 3 Spatial distributions of multi-year average annual precipitation (a), annual average air temperature (b), 
soil moisture (c), and land use and land cover change (d) in the HRB 


In terms of drought propagation probability, both seasonal sensitivity and drought level sensitivity 
(i.e., the probabilities of agricultural drought at different levels triggered by meteorological 
drought at different levels) were investigated. 

The framework of this paper is presented in Figure 4. In particular, we focus on the following 
three issues: (1) the spatial characteristics of meteorological drought and agricultural drought 
across the basin from 1981 to 2020; (2) the seasonal sensitivity and dynamic characteristics of 
drought propagation time; and (3) both the seasonal sensitivity and drought level sensitivity of 
drought propagation probability. 
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Fig. 4 Framework of this study. SPEI, standardized precipitation evapotranspiration index; SWDI, soil water 
deficit index. SPEI-1, SPEI-3, SPEI-6, and SPEI-12 mean SPEI at the 1-, 3-, 6-, 12-month time scales, 
respectively. 


Meteorological drought and agricultural drought were characterized by SPEI and soil water 
deficit index (SWDI), respectively. As SPEI has the characteristics of multiple time scales, we 
used SPEI at the 1-, 3-, 6-, and 12-month time scales (denoted as SPEI-1, SPEI-3, SPEI-6, and 
SPEI-12, respectively) to investigate the spatial patterns of meteorological drought. The seasonal 
sensitivity and dynamic characteristics of drought propagation time were evaluated using the 
correlation coefficient method, sliding windows, and Mann-Kendall test. Propagation probability 
and its sensitivities to seasons and drought levels were quantified by the Copula-based interval 
conditional probability model. All data used here were from GLDAS-2.0 and GLDAS-2.1. 
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Unlike the regular four seasons in one year (spring, summer, autumn, and winter), we only 
consider the following two seasons in this study due to the high frequency of irrigation in the 
study area: (1) the irrigation season from April to September; and (2) the non-irrigation season 
from October to March of the following year. 


3.2 Drought indices 


3.2.1 Standardized precipitation evapotranspiration index (SPEI) 

SPEI, proposed by Vicente-Serrano et al. (2010), is a drought index based on the standardized 
precipitation index (SPI) by introducing potential evapotranspiration. The Thornthwaite approach, 
which requires only monthly precipitation and air temperature data, is suggested to calculate the 
potential evapotranspiration (Li et al., 2016). SPEI is obtained after normalization by Equations 1 
and 2: 


Co +¢,@+c,0" 
+ | 
1+d,o+d,o0 +d,0° 


@=./-2InF(D) (0<F(D)< 0.5), (2) 
where œ is the function of F(D), and F(D) is the probability distribution function of D, which is 
the difference series between precipitation and potential evapotranspiration; and co, ci, c2, di, d2, 
and d; are constants, with the values of 2.515517, 0802853, 0.010328, 1.432788, 0.189269, and 
0.001308, respectively. Further details of the calculation process can be found in Vicente-Serrano 
et al. (2010). Table 2 shows the classification of meteorological drought categories according to 
SPEI values. 


SPEI =a (1) 


Table 2 Classification of meteorological drought categories based on standardized precipitation evapotranspiration 
index (SPEI) 


Level Meteorological drought category SPEI values 
0 Non-drought SPEI> -0.5 
1 Mild drought —1.0<SPEI< —0.5 
2 Moderate drought —1.5<SPEI< -1.0 
3 Severe drought —2.0<SPEI< -1.5 
4 Extreme drought SPEI< -2.0 


3.2.2 Soil water deficit index (SWDI) 


Unlike many drought indices based on meteorological variables or water balance equations, 
SWDI is based on a previous water deficit index (Martinez-Fernandez et al., 2015). It considers 
the relationship between plant physiological state and soil moisture, and is more suitable for 
monitoring agricultural drought compared to other indices (Zhu et al., 2019; Wu et al., 2020). 
SWDI is calculated as follows: 


0-6 
swoi =| Fe )x1o, (3) 


Pxwe = Orc — Awe > (4) 
where @ is the soil moisture (g/cm*); and Orc, Owe and awc represent the field capacity (g/cm*), 
wilting point (g/cm*), and available water capacity (g/cm°), respectively. The long-term minimum 
and maximum values based on the soil moisture time series are estimators of Orc and Owp, 
respectively (Bai et al., 2018). 

Table 3 lists the agricultural drought categories based on the values of SWDI 
(Martinez-Fernandez et al., 2015). 
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Table 3 Classification of agricultural drought categories based on soil water deficit index (SWDI) 


Level Agricultural drought category SWDI values 
0 Non-drought SWDL>0 
1 Mild drought —2<SWDISO 
2 Moderate drought —5<SWDIk< -2 
3 Severe drought —10<SWDI< -5 
4 Extreme drought SWDIs -10 


3.3 Determination of drought propagation time 


3.3.1 Pearson correlation coefficient 


Propagation time can be defined as the time from the onset of meteorological drought to the 
beginning of agricultural drought, and can be simplified into a mathematical linkage between 
SPEI and SWDI in different accumulation periods. The time scale for meteorological drought is n 
months (n=1, 2, ..., 12), while the time scale for agricultural drought is 1 month. The maximum 
Pearson correlation coefficient between SPEI-n and SWDI-1 represents the connection between 
meteorological drought and agricultural drought, and the accumulation period n is defined as the 
drought propagation time (Xu et al., 2021). 

3.3.2 Mann-Kendall test and sliding windows 


The Mann-Kendall test is a nonparametric statistical test used to quantitatively detect any 
significant trends in long term series. It does not require samples to follow a certain distribution 
and is less affected by outliers; thus, the Mann-Kendall test has been widely used in trend and 
mutation testing in the fields of hydrology and meteorology (e.g., Sun, 2021). The Sen's slope 
estimator is often used together with the Mann-Kendall test to quantify the magnitude of change. 
In particular, the Mann-Kendall test is applied to test the direction of the trend, while the Sen's 
slope estimator is employed to determine the magnitude of the trend (Bazrafshan, 2017; Moral et 
al., 2017). 

The sliding window algorithm includes two important elements: the sliding distance and 
window length. By assigning these two element values, the new series is taken in a certain 
direction according to a fixed sliding distance and length (Shao et al., 2019). In this study, the 
sliding distance was set as 1 a and the sliding window was set as 20 a. Based on these settings, we 
divided the index series from 1981 to 2020 into 21 new series, and the dynamic changes of 
drought propagation time were then determined by calculating the values of the Sen's slope 
estimator for these 21 series. 


3.4 Copula function and conditional probability 


There are two commonly used Copula functions: Archimedean Copula and Elliptical Copula 
(Sun, 2021). The Archimedean Copula is highly representative and easy to construct, while the 
Elliptical Copula can establish non-normal relationships and better describe extreme events. The 
Copula function is calculated as follows: 

C(F (u), F, 0) =F (F(u), PO), (5) 
where C(Fı(u), F2(v)) is the Copula function combined with two random variables u and v; F 
denotes the marginal distribution function; and F; and F> represent the marginal distribution 
functions of the two random variables, respectively. 

Five Copula functions (including four Archimedean Copula functions and one Elliptical Copula 
function) were selected to establish the joint distribution of the meteorological drought and 
agricultural drought indices. Table 4 details their formulas and parameters. The Akaike 
information criterion (AIC) and Bayesian information criterion (BIC) were used to evaluate the 
goodness-of-fit of the functions. The Copula function with the smallest AIC or BIC value is 
considered as the best-fit function for the SPEI and SWDI series. Before establishing the Copula 
function, five alternative marginal distribution functions (i.e., normal, logistic, Weibull, general 
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extreme value (GEV), and Pearson-III (P3)) were chosen to fit the SPEI and SWDI series. The 
Kolmogorov-Smirnov (K-S) test was used to evaluate the performance of each marginal 
distribution function. 


Table 4 Cumulative density functions and parameter ranges of Copula functions 


Copula function type C(u, v) r value range 
Clayton (ur tv" i r>0 
Gumbel exp(-((-lnu)" + (-Invy’ yr) r21 
Archimedean Copula m a _ De” = ») = 
7 e” -l1 
a 1-(i-u)' +(1-v)! ~a-wid-wy Fal 


a a 2 2 _k+2 
+ ts CH) pRO) S° — 2rst +t" ast 
Elliptical Copula t f ‘ f i : 
a 


1 
exp[l+ 
> 2my1-r° k(l-r°) 


Note: u and v are two variables; r is the parameter; f,' is the inverse function of the univariate ¢ distribution function with a 
degree of freedom of k; and s and ¢ are the univariate ¢ distribution functions of variables u and v. 


An interval conditional probability model can characterize the relationship between 
conditionally restricted random variables (Sattar et al., 2019). The interval conditional probability 
models of SPEI and SWDI were used to explore the probability of different levels of agricultural 
drought caused by different levels of meteorological drought. The corresponding formulas are 
described as: 


P(v <V < v,,u SU <u) 
P(u SU <u) 
a CE (u), Fa 0) — Ciu), Fo (%)) -Ci 2), )) + CH U2), Foa) 
F, (uz) — Fu) 

where P is the conditional probability value; U and V represent the SPEI and SWDI series, 
respectively; vı and uı are the lower limits of variables v and u, respectively; v2 and u2 are the 
upper limits of variables v and u, respectively; Fı and Fz represent the marginal distribution 
functions of variables u and v, respectively; and C(Fi(w1), F2(v1)), C(Fi(u1), F2(v2)), CLF (v2), 
F(v1)), and C(Fı(u2), F2(v2)) are the corresponding Copula functions. 


Pv, <V <v, |u <U <u,)= 


(6) 


4 Results 


4.1 Meteorological drought and agricultural drought over the study area 


Figure 5a—d show the spatial distribution of multi-year average meteorological drought based on 
SPEI over the HRB at different time scales from 1981 to 2020. The study area was generally 
dominated by the non-drought category at the multi-year average scale, although SPEI values 
differed with time scales. Larger SPEI values were observed based on SPEI-1 (mainly in green 
and yellow colors) with only 1% negative values, while smaller values were observed based on 
SPEI-3, SPEI-6, and SPEI-12 (mainly in orange and red colors), exhibiting much more negative 
values (37%, 31%, and 25% of the total values, respectively). 

Figure 5e exhibits the spatial distribution of agricultural drought characterized by SWDI. The 
blue and green colors indicate non-drought with SWDI values greater than zero, while the yellow 
and red colors indicate drought with SWDI values lower than zero. With the exception of the river 
system and its surroundings in the midstream and downstream regions, the majority of the study 
area was dominated by orange, light green, and red colors, indicating the presence of agricultural 
drought in these regions. Extreme agricultural drought was detected in the eastern downstream 
region (red color in Fig. 5e), which is probably related to the low soil moisture in this region (Fig. 
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3c). Extreme agricultural drought was also detected in some upstream regions (orange color in 
Fig. 5e), although the intensity was not as severe as those in the eastern downstream region. The 
river system and its surrounding areas (blue color in Fig. 5e) were categorized as non-drought, 
which may be related to the high soil moisture. 
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Fig. 5 Spatial distributions of SPEI-1 (a), SPEI-3 (b), SPEI-6 (c), SPEI-12 (d), and SWDI (e) in the HRB 
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4.2 Drought propagation time 


4.2.1 Drought propagation time and its seasonal sensitivity 

Figure 6 shows the spatial distribution of drought propagation time for the year-round, irrigation 
season, and non-irrigation season in the HRB. The redder the color tone, the longer the drought 
propagation time. Drought propagation time was observed to vary in different regions of the 
basin. Taking the case of Figure 6a as an example, the midstream and southern downstream 
regions experienced a short drought propagation time (less than 2 months), while other regions 
experienced a longer drought propagation time for the year-round. This is particularly true for the 
upstream and northern downstream regions, with a drought propagation time of more than 11 
months. The average drought propagation time in the whole basin reached about 8 months, with 
the upstream region having the longest propagation time (an average of 12 months), followed by 
the midstream region (8 months), and the downstream region (6 months) (Fig. 6a). 

Figure 6b and c reveals that drought propagation time exhibited seasonal characteristics. More 
specifically, drought propagated much faster and took less than 2 months to transform from 
meteorological drought to agricultural drought during the irrigation season in the midstream and 
downstream regions. In contrast, it lasted for more than 7 months during the non-irrigation 
season. On average, drought propagation time was shorter during the irrigation season, with an 
average time of 6 months in the whole basin, and 9, 5, and 3 months in the upstream, midstream, 
and downstream regions, respectively (Fig. 6b). However, during the non-irrigation season, 
drought propagation time became longer, with an average time of 10 months in the whole basin, 
11 months in both the upstream and midstream regions, and 9 months in the downstream region 
(Fig. 6c). 
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Fig. 6 Spatial distributions of propagation time from meteorological drought to agricultural drought for the 
year-round (a), and during the irrigation season (b) and non-irrigation season (c) in the HRB 


4.2.2 Dynamic characteristics of drought propagation time 


To explore the dynamic characteristics of drought propagation time, we set a sliding window of 
20 a and divided the study period (1981-2020) into 21 time series. Following this, the trend of 
drought propagation time in each grid was detected using the Sen's slope estimator and 
Mann-Kendall test. Figure 7 depicts the spatial distribution of the temporal trends for drought 
propagation time. All tests exceeded the 95% confidence level, indicating that the trends are 
significant at the 0.05 significance level. The green color indicates the Sen's slope estimator 
values greater than zero (upward trend), while the red color means the values less than zero 
(downward trend). 

From Figure 7a it can be seen that at the year-round time scale, the upstream and northern 
downstream regions (red color; accounting for 56% of the basin) exhibited downward trends of 
drought propagation time. This indicates a shortening of drought propagation time in these 
regions. The midstream and southern downstream regions (green color; accounting for 44% of the 
basin) showed upward trends, revealing that drought propagation time was longer in these 
regions. 
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Fig. 7 Spatial distributions of the Sen's slope estimator values (P<0.05) of propagation time from 
meteorological drought to agricultural drought for the year-round (a), and during the irrigation season (b) and 
non-irrigation season (c) in the HRB 


During the irrigation season, more regions (78% of the basin) presented upward trends (green 
color in Fig. 7b). However, opposite trends were observed during the non-irrigation season, and 
nearly 98% of the basin exhibited downward trends (red color in Fig. 7c), indicating that the 
propagation time from meteorological drought to agricultural drought in most regions of the basin 
became shorter. More specifically, with the progress of time, agricultural drought triggered by 
meteorological drought occurred more quickly during the non-irrigation season than during the 
irrigation season. 

4.3 Drought propagation probability 

4.3.1 Selection of optimal marginal distribution and copula functions 

Drought propagation probability was calculated using the Copula-based interval conditional 
probability model. The normal, logistic, Weibull, GEV, and P3 marginal distribution functions 
were used to fit the SPEI and SWDI series. Figure 8 shows the percentage of each marginal 
distribution function in all grids. GEV ranked the first for the SPEI-1, SPEI-6, and SPEI-12 series 
(with the percentage values of 64%, 48%, and 54%, respectively) and the second for the SPEI-3 
series (with a percentage value of 28%). For SWDI, P3 ranked the first. In order to maintain 
consistency in selecting an optimal marginal distribution function for each type of drought index, 
we selected GEV and P3 as the marginal distributions for the SPEI (including SPEI-1, SPEI-3, 
SPEI-6, and SPEI-12) and SWDI series, respectively. 
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Fig. 8 Percentages of different marginal distribution functions for the SPEI and SWDI series. GEV, General 
extreme value; P3, Pearson-MI. 


To ensure the reliability of the selected marginal distribution functions, we also performed the 
K-S test for each grid. For the SPEI series, more than 94% of the grids can be fitted by GEV, and 
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for the SWDI series, 82% of the grids can be fitted by P3, according to the K-S test results shown 
in Table 5. Thus, GEV and P3 were finally determined as the optimal marginal distribution 
functions for the SPEI and SWDI series, respectively. 


Table 5 Percentages of grids that passed the Kolmogorov-Smirnov (K-S) test for the optimal marginal 
distribution function 


Drought index Optimal marginal distribution function Percentage of grids that passed the K-S test (%) 
SPEI-1 GEV 99.95 
SPEI-3 GEV 99.98 
SPEI-6 GEV 98.13 
SPEI-12 GEV 94.07 
SWDI P3 81.98 


Note: GEV, General extreme value; P3, Pearson-III. 


The Copula functions were fitted using R (v.4.1.1) (R Development Core Team, 2022), and the 
function with the largest percentage among all grids was selected as the optimal one. Figures 9 
and 10 show the percentage of each Copula function under the AIC and BIC. The Clayton Copula 
function (light pink in Fig. 9 and green color in Fig. 10) accounted for the highest percentage for 
the combinations of SPEI-1-SWDI, SPEI-3-SWDI, and SPEI-6-SWDI, and the second-highest 
percentage for SPEI-12-SWDI (slightly inferior to the Gumbel Copula function). Thus, the 
Clayton Copula function was finally chosen as the optimal Copula function for the SPEI and 
SWDI series. 
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Fig. 9 Percentages of Copula functions for the combinations of SPEI and SWDI series among all grids in the 
HRB under the Akaike information criterion 


60 


Clayton =Gumbel = Frank m Joe 


50 


40 


‘LLU 


SPEI-1-SWDI  SPEI-3-SWDI  SPEI-6-SWDI_SPEI-12-SWDI 
Combination of drought index 


Percentage (%) 


Fig. 10 Percentages of Copula functions for the combinations of SPEI and SWDI series among all grids in the 
HRB under the Bayesian information criterion 
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4.3.2 Drought propagation probability and its drought level sensitivity 
Drought propagation probability was subsequently calculated using Equations 5 and 6 and the 
Clayton Copula function. In addition, SPEI-6 was selected to calculate drought propagation 
probability. This is because the average propagation time from meteorological drought to 
agricultural drought in the HRB was estimated to be 8 months, which is close to 6 months. In 
addition, such a medium time scale is suitable for characterizing the propagation of water deficit 
signals from the atmosphere to the soil (Xu et al., 2021; Yang et al., 2021; Zhang et al., 2021b). 
Figure 11 shows the propagation probability of agricultural drought at different levels induced 
by varying levels of meteorological drought. Drought propagation probability ranged from 0.00 to 
1.00, and the redder the color tone, the greater the drought propagation probability value. 
Generally speaking, agricultural drought caused by meteorological drought exhibited varying 
spatial trends. 
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Fig. 11 Spatial distributions of the propagation probability values of different agricultural drought levels (i.e., 
mild, moderate, severe, and extreme drought levels from the top to bottom) under various meteorological drought 
levels (i.e., mild, moderate, severe, and extreme drought levels from the left to right) at the year-round scale in the 
HRB 


For the upstream region, when meteorological drought occurred, the propagation probability of 
severe agricultural drought was about 0.60 (mainly in blue color in Fig. 11i and j), and that of 
extreme agricultural drought was 0.40—0.85 (green and purple colors in Fig. 11m and n), which 
was slightly greater than those of mild and moderate agricultural droughts, with the propagation 
probability values of 0.20 and 0.40 (orange color in Fig. lla and b for mild agricultural drought, 
and green color in Fig. lle and f for moderate agricultural drought). For the eastern upstream 
region, extreme agricultural drought was more likely to be triggered (Fig. 11m—p), with 
propagation probability values higher than 0.80. The propagation probability of agricultural 
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drought at different levels increased with the meteorological drought levels. For example, when 
the level of meteorological drought increased from moderate to severe, the propagation 
probability of agricultural drought increased from 0.40 (moderate level) to 0.60 (severe level). 

For the midstream region, the propagation probability of agricultural drought at different levels 
also increased with the meteorological drought levels. For example, the propagation probability of 
severe agricultural drought increased from 0.80 to over 0.90 when the meteorological drought was 
from moderate level to severe level (mainly from purple to red color in Fig. 11j and k). For the 
upstream-midstream junction region, severe and moderate agricultural drought was more likely to 
occur (Fig. 11e—l), with propagation probability values higher than 0.80. 

The downstream region was mainly dominated by mild, moderate, and extreme agricultural 
droughts induced by meteorological drought. Extreme drought was most likely to occur in the east 
of the river system (Fig. 11m-p), while moderate agricultural drought was more probable in the 
west of the river system, with propagation probability values greater than 0.90 (Fig. 11e—h). For 
the southern downstream region (near the midstream region), the propagation probability of mild 
and moderate agricultural droughts increased with the meteorological drought levels. For 
example, the propagation probability of moderate agricultural drought increased from 0.40 to 0.80 
and subsequently exceed 0.90 when the meteorological drought increased from the mild level to 
moderate level and then to the extreme level (Fig. 1le—h). Note that the propagation probability 
values near the river system were extremely close to zero. This is because soil moisture near the 
mainstream region was always high, implying that the propagation probability of agricultural 
drought tends to be low. 

In short, with the exceptions of regions close to the river, different agricultural drought levels in 
the HRB were sensitive to different meteorological drought levels, and such sensitivity exhibited 
noticeable spatial differences. 

4.3.3 Seasonal sensitivity of drought propagation probability 

Figures 12 and 13 present the spatial distributions of the seasonal sensitivity of drought 
propagation probability over the basin. Meteorological drought was more likely to trigger 
agricultural drought at low levels during the irrigation season and at severe level during the 
non-irrigation season. Such seasonal sensitivity of drought propagation probability also had 
noticeable spatial differences. 

In particular, for the western downstream region, only mild and moderate agricultural droughts 
were triggered with the propagation probability values of 0.20-0.40 and >0.90 during the 
irrigation season, respectively (Fig. 12a—h), while moderate and severe agricultural droughts were 
triggered during the non-irrigation season, with both propagation probability values higher than 
0.85 (Fig. 13a—h). For the southern downstream region, only mild agricultural drought was 
triggered, with propagation probability values over 0.90 during the irrigation season (Fig. 12a—d), 
while mild and moderate agricultural droughts were triggered during the non-irrigation season, 
with propagation probability values greater than 0.85 (Fig. 13a—d). For the eastern downstream 
region, extreme agricultural drought was most likely to occur with propagation probability 
exceeding 0.90 regardless of the irrigation and non-irrigation seasons. 

For the midstream region, mild, moderate, severe, and extreme agricultural droughts (probability 
values ranging from 0.20 to higher than 0.90) occurred during both the irrigation and 
non-irrigation seasons. However, the area of severe agricultural drought triggered during the 
non-irrigation season was larger than that during the irrigation season (Figs. 12i-l and 13i-l). As 
for the upstream region, minimal changes were detected for the propagation probability of 
agricultural drought with varying levels during both the irrigation and non-irrigation seasons 
(Figs. 12i-] and 13i-1). 


5 Discussion 


It is well known that there is a water and energy link between meteorological drought and 
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Fig. 12 Spatial distributions of the propagation probability values of different agricultural drought levels (i.e., 
mild, moderate, severe, and extreme drought levels from the top to bottom) under various meteorological drought 
levels (i.e., mild, moderate, severe, and extreme drought levels from the left to right) during the irrigation season 
in the HRB 


agricultural drought, with a phase difference in time. The propagation between them is 
inseparable from the water cycle (Xu et al., 2021; Li et al., 2022). In this study, we analyzed the 
propagation time and propagation probability from meteorological drought to agricultural 
drought, as well as their sensitivities to seasons and drought levels in the HRB. The drought 
propagation characteristics and corresponding sensitivities exhibited great spatial differences, 
which may be related to the local climate, topography, vegetation cover, etc. (Xu et al., 2019; 
Zhou et al., 2021a; Zhu et al., 2021; Zhang et al., 2022b). 


5.1 Drought propagation time 


The upstream region in the HRB was observed to experience a longer propagation time from 
meteorological drought to agricultural drought, which is probably related to the land use and land 
cover type. As shown in Figure 3d, the upstream region is mainly covered by forest vegetation, 
which has a good water holding capacity. In the short term, changes in precipitation have a 
limited effect on soil moisture. This can result in a slow response of agricultural drought to short 
term precipitation changes, and thus a longer propagation time from meteorological drought to 
agricultural drought in the upstream region. 

In the midstream and southern downstream regions, drought propagation time was only 1-2 
months. This may be related to the climate and geographical conditions of these regions. The 
midstream of the basin is an agricultural planting region that requires a lot of irrigation water to 
maintain soil moisture and ensure crop growth. The southern downstream region is dominated by 
marshland, saline land, and desert, and soil moisture is relatively low, making it particularly 
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Fig. 13 Spatial distributions of the propagation probability values of different agricultural drought levels (mild, 
moderate, severe, and extreme drought levels from the top to bottom) under various meteorological drought levels 
(i.e., mild, moderate, severe, and extreme drought levels from the left to right) during the non-irrigation season in 
the HRB 


important for precipitation to replenish soil moisture. Thus, these regions are more sensitive to 
changes in climatic conditions. If precipitation continues to be scarce (without considering the 
artificial diversion irrigation), the soil moisture will be reduced, thus quickening the response of 
agricultural drought (characterized by soil moisture) to meteorological drought. The northern 
downstream region contains an inland lake and grassland, where soil moisture is relatively high. 
Furthermore, the average annual precipitation in this region is minimal (less than 50.0 mm) (Ding 
et al., 2009; Yang et al., 2022). These may lead to a weak sensitivity of this region to the changes 
in climatic conditions and thus a slow response of agricultural drought to meteorological drought. 
Other factors, including wind, radiation, and altitude, could also affect soil moisture availability 
by influencing evapotranspiration and runoff processes (Van den Hoof and Lambert, 2016). 

In the HRB, the irrigation season tends to experience a shorter propagation time from 
meteorological drought to agricultural drought than the non-irrigation season. This indicates that 
agricultural drought is more sensitive to meteorological drought during the irrigation season. 
More than 70% of precipitation is concentrated in June—-September (during the irrigation season), 
and if the lack of precipitation continues to occur, the soil moisture will decrease accordingly. In 
addition, the period from June to August is the growing season and vegetation growth requires 
more water from the soil (Qian et al., 2016). Agricultural drought is more likely to be triggered by 
the lack of precipitation and high-water demand for vegetation growth. Thus, there is a relatively 
shorter drought propagation time during the irrigation season. 

Figure 6b and c reveals that, irrespective of the irrigation or non-irrigation season, the upstream 
region always experienced a longer drought propagation time compared to other regions of the 
basin. This may be because the upstream region is dominated by forest, and the roots of the 
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vegetation are deep and well developed, thus leading to higher resistance to the reduction of 
surface soil moisture compared to crops, grassland, and marshland covered in the midstream and 
downstream regions. Moreover, drought propagation time during the irrigation season in the 
upstream region was slightly shorter than that during the non-irrigation season, which may be due 
to the presence of snow and ice cover in the basin (Chen et al., 2018; Li et al., 2021). As 
temperatures increase, snow and permafrost continue to melt, replenishing soil moisture (LYU et 
al., 2020) and thus shortening drought propagation time during the irrigation season. 

Although few studies have investigated the seasonal sensitivity of propagation time from 
meteorological drought to agricultural drought over the HRB, the relevant cases in other basins 
can provide reference for us (e.g., Huang et al., 2015; Bhardwaj et al., 2020; Li et al., 2022). For 
example, Huang et al. (2015) concluded that the propagation time from meteorological drought to 
agricultural drought in the Weihe River Basin is shorter in summer, and longer in autumn. Liang 
et al. (2021) estimated that the propagation time from meteorological drought to agricultural 
drought in the Jinta River Basin of Northwest China was concentrated in 2 months during 
summer, and 6—12 months in autumn and winter. Li et al. (2022) found that the propagation time 
from meteorological drought to agricultural drought was shorter (1-3 months) in summer and 
autumn while longer (5-12 months) in spring and winter over an arid region of Northeast Asia. 
The previous literature shows that the propagation time from meteorological drought to 
agricultural drought in summer is relatively shorter, which is similar to our findings in the HRB 
during the irrigation season. 


5.2 Drought propagation probability 


The propagation probability of agricultural drought at different levels increased with the 
meteorological drought levels in the upstream and midstream regions of the HRB. This may be 
related to the vegetation cover and evapotranspiration. When mild meteorological drought occurs, 
the excellent water retention capacity of the vegetation can withstand the effects of 
meteorological drought. However, an extremely reduction in precipitation can induce a 
water-stressed state of the vegetation, consuming more water through evapotranspiration during 
droughts, which may exacerbate droughts and induce a higher propagation probability of 
agricultural drought (Zhu et al., 2021). The propagation probability of agricultural drought 
triggered by meteorological drought in the midstream region was always higher than that in the 
upstream region, which is probably due to the different vegetation cover types. The forest in the 
upstream region can be more drought tolerant because the well-developed root systems can 
absorb water from deeper soils, while grasses and crops in the midstream region have shorter and 
undeveloped root systems, which require more shallow soil water to support the growth of them 
(Zhou et al., 2021a). For the downstream region, the levels of agricultural drought varied at the 
spatial scale, which was inferred to be largely related to soil moisture. The low soil moisture in 
the western and eastern downstream region increased the propagation probability of moderate and 
extreme agricultural droughts, while the high soil moisture near the river system resulted in 
agricultural drought less likely to be triggered by meteorological drought. 

At the seasonal scale, lesser agricultural droughts were more likely to occur during the 
irrigation season in the midstream and downstream regions, while severer agricultural droughts 
mostly occurred during the non-irrigation season. This implies that irrigation-induced soil 
moisture increases could delay and mitigate the effects of meteorological drought on agricultural 
drought. Irrigation affects the water and energy balance near the ground, making the atmosphere 
wetter near the surface (Zhu et al., 2021). Therefore, it is speculated that the same level of 
meteorological drought is more likely to trigger a lower level of agricultural drought during the 
irrigation season. This result is in agreement with the findings of Fang et al. (2019). In contrast, 
no obvious differences of propagation probability were observed in the upstream region for both 
the irrigation and non-irrigation seasons (Figs. 12 and 13) since soil moisture in the upstream 
region was not affected by irrigation. 
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6 Conclusions 


This study investigated the propagation time and propagation probability from meteorological 
drought to agricultural drought in the HRB from 1981 to 2020, together with their spatial 
variations and sensitivities to seasons and drought levels. The correlation coefficient method and 
Copula-based interval conditional probability model were employed for the analysis. The 
following conclusions can be drawn from this study. 

(1) Larger SPEI values were observed based on SPEI-1, while smaller SPEI values were 
observed based on SPEI-3, SPEI-6, and SPEI-12. Spatially, smaller SPEI values occurred in the 
midstream and northern downstream regions. Agricultural drought characterized by SWDI 
exhibited the most serious level in the eastern downstream region, followed by the upstream 
region. 

(2) The average drought propagation time was estimated to be 8 months in the whole basin at 
the year-round scale, 6 months during the irrigation season, and 10 months during the 
non-irrigation season. Drought propagation time had obvious spatial differences in the HRB. By 
analyzing its dynamic characteristics, we found an increasing trend in 78% of the basin during the 
irrigation season, and a decreasing trend in 98% of the basin during the non-irrigation season. 

(3) Propagation probability was sensitive to seasons and drought levels, and such sensitivity 
had noticeable spatial differences in the whole basin. In the upstream, midstream, and southern 
downstream regions, the propagation probability of agricultural drought at different levels 
increased with the meteorological drought levels. For the western and eastern downstream region, 
moderate and extreme agricultural droughts were detected to be more easily triggered by 
meteorological drought, with propagation probability values higher than 0.90. In terms of 
seasonal sensitivity, lesser agricultural droughts were more likely to be triggered by 
meteorological drought during the irrigation season, while severer agricultural droughts were 
occurred mostly during the non-irrigation season. 

The results can help to improve our understanding of drought propagation among different 
drought types, and can act as a useful reference for monitoring meteorological drought and 
agricultural drought and the rational use of water resources in the HRB. However, this study also 
has some limitations. For example, the influencing factors of drought propagation characteristics 
were only explained through qualitative analysis. Quantitative analysis of climate, land use, soil 
moisture, and irrigation on drought propagation will be the focus of the future work. 
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